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The Purcell effect - the modification of the spontaneous emission rate in presence of resonant 
cavities or other resonant objects - is a fundamental effect of quantum electrodynamics. However, 
a change of the emission rate caused by environment different from free space has a classical coun¬ 
terpart. Not only quantum emitters, but any small antenna tuned to the resonance is an oscillator 
with radiative losses, and the influence of the environment on its radiation can be understood and 
measured in terms of the antenna radiation resistance. We present a general approach which is 
applicable to measurements of the Purcell factor for radio antennas and to calculations of these 
factors for quantum emitters. Our methodology is suitable for calculation and measurement of both 
electric and magnetic Purcell factors, it is versatile and applies to various frequency ranges. The 
approach is illustrated by a general equivalent scheme and allows the Purcell factor to be expressed 
through the continious radiation of a small antenna in presence of the environment. 
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I. INTRODUCTION 

The Purcell effect is defined as a modification of the 
spontaneous emission lifetime of a quantum source in¬ 
duced by the interaction with its environment [1-10]. 
This modification is significant if the environment is a 
resonator tuned to the emission frequency. This effect 
has been first pointed out by E.M. Purcell [1] in 1946 in 
the context of nuclear magnetic resonance experiments. 
At present, this effect is widely used in microcavity light- 
emitting devices [11-13], in single-molecule optical mi¬ 
croscopy [14-17], in microscopy of single NV centers in 


nanodiamonds [18], of Eu 3+ -doped nanocrystals [8] and 
for visualization of biological processes with participa¬ 
tion of large molecules, such as DNA [19]. An overview 
of nanosensing applications of the Purcell effect is pre¬ 
sented in [20]. 

We start from a brief overview of several equivalent 
definitions of the Purcell factor - the value which de¬ 
scribes this effect quantitatively - and the existing ap¬ 
proaches to its theoretical and experimental evaluation. 
First, let us point out that the Purcell effect is based 
on the quantum electrodynamics concept of weak cou¬ 
pling of an emitter and a resonating object (nanoan¬ 
tenna [20], or optical cavity [11-13]). The weak and 
strong coupling regimes [21, 22] can be distinguished by 
comparing the so-called emitter-field coupling constant 
X = [|d| 2 o;o/(2feoV)] 1 / 2 with the decay rate of the pho¬ 
ton in a cavity 7 and the nonradiative decay rate of the 
excited state qais- Here, to 0 and d = e(2|r|l) are the fre¬ 
quency of the excited-to-ground state transition (2 —> 1) 
and its dipole moment (matrix element), respectively, e 
is the electron charge, V is the effective volume of the 
resonator mode, e 0 is the vacuum permittivity. We use 
the SI units, the result in the CGS units can be obtained 
by replacing e 0 by 1/(47t). 

In the weak-coupling regime, when x "C 7,7dis: the 
hybridization of the quantum emitter and the resonator 
eigenstates is weak. Therefore the frequency w 0 of the 
spontaneous emission is not modified by the resonator, 
and the interaction only leads to a modification of the 
decay rate. The dipole moment of the optical transition 
d and its classical dipole moment di keep unperturbed 
and di = 2d (see for example Ref. [23], pp.250-251). The 
ratio of the decay rate 7 in the vicinity of the resonator 
to the decay rate of the same emitter in free space 70 can 
be written as [21]: 
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Figure 1. The classical (a) and quantum (b) realizations of the Purcell effect. 


where 70 = wd|di| 2 / (12TT£ 0 hc 3 ) = wjj|d| 2 /(37reo?ic 3 ) [21], 
q = Lu/c is wavenumber in free space, E s (rd) is the mag¬ 
nitude of an electric field of the quantum source dipole 
di oscillating at the frequency uiq, scattered from an in¬ 
homogeneous environment and evaluated at the source 
origin r^. The quantity P is called the Purcell factor. 
According to Eq.(l), the magnitude of the Purcell fac¬ 
tor does not depend on the magnitude of the transition 
dipole moment d, because the scattered field value is di¬ 
rectly proportional to the dipole moment. For the pur¬ 
pose of this paper it is important to note, that equa¬ 
tion (1) can be applied not only to the cavities or nanoan¬ 
tennas but to an arbitrary electromagnetic environment 
of the emitter different from free space [20]. Moreover, 
the concept of Purcell’s factor can be extended to opti¬ 
cal emitters which cannot be modeled as a point electric 
dipole [24, 25]. The Purcell factor can be also understood 
in terms of the local density of photonic states modified 
by the presence of the object [26]. 

The above expression for 70 does not take into account 
the non-radiative decay (it is assumed that 7di s < 7o) 
and results from the standard formula for the power ra¬ 
diated by a Hertzian dipole di at frequency u> q: 
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and the Purcell factor can be expressed also as 


I? _ T _ ° ra d + Pnonrad _ p , / c \ 

F 0 — -^rad 1 -^nonrad- 

70 -M3,rad 


Here, P r ad is the power radiated in the far zone (enhanced 
by the environment) and P n0 nrad is the power dissipated 
in the environment. 

If the electromagnetic environment is lossless, the last 
term vanishes and the Purcell factor describes the change 
of the total radiated power P ra( j at the frequency of the 
emitter: 


p = 2L — ^r ad 

70 To,rad 


( 6 ) 


where the index 0 still means the corresponding value for 
the same emitter in free space. If the emitter is located 
in a lossy medium (perhaps inhomogeneous) the Purcell 
factor Eq. (1) has two contributions: that corresponding 
to the far-held emission and that corresponding to the 
Joule losses in the environment [27]. When this environ¬ 
ment can be described by position-dependent dielectric 
constant e(r') the Joule loss contribution into Purcell’s 
factor can be presented as [28]: 

Pnonrad = ^^2 j d V Im^r')] |E(r')| 2 , (7) 


namely, 
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III.rad 
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is the ratio of Po, ra d to the photon energy. In the weak 
coupling regime, the environment modifies only the radi¬ 
ated (far-zone) power and the dissipation of power in the 
volume outside of the emitter. Thus, the decay factor 
modified by the environment can be written as 


where E(r') is the total held produced by the dipole d\ 
at the point r' which is integrated over the surrounding 
space. 

Here, one may introduce the radiation efficiency of the 
quantum source £ in the same way as it is done in the 
antenna theory [29]: £ = P ra d/P = (P - P n 0 nrad)/P- 
The total quantum yield of the emitter Q is determined 
by the competition between the far-held radiation, the 
Joule losses, and the internal non-radiative losses of the 
emitter 7^: 


7 = 


Prad T Pnonrad 


hujQ 
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Note, that in formula (6) we assumed the decay rate 
of the emitter in free space to be equal 70, i.e. neglected 
the non-radiative losses inside the emitter. This can be 
a realistic approximation for many quantum dots and 
fluorescent dye molecules (e.g. in [30] 7 d i s and 70 were 
separately measured for nanocrystal quantum dots and 
it was shown that 7di s < 70). 

For quantum emitters the total Purcell factor is mea¬ 
sured either directly by evaluating the speedup of the 
time-resolved photoluminescence [9] or indirectly, for ex¬ 
ample, using the Raman spectroscopy [31]. High values of 
F can be achieved with nanoantennas - resonant devices 
that effectively convert the near field of quantum sources 
to propagating optical radiation [32, 33]. This transfor¬ 
mation is carried out by means of impedance matching 
between the quantum source and the nanoantenna [34— 
37]. Another possibility to attain large values of the Pur¬ 
cell factor Eq. (1) is provided by hyperbolic metamateri¬ 
als (see the review [9]). 

In this paper, we look at the Purcell effect in the 
broader context and study it within the classical frame¬ 
work. While the question whether the spontaneous decay 
itself is a truly quantum [38] or not a purely quantum 
phenomenon [39, 40] is still under debate, the modifica¬ 
tion of the spontaneous decay in a medium can be def¬ 
initely considered classically. Indeed, the electric field 
E d (r d ) of the emitter (’’antenna” in the Fig.1(a) or 
’’quantum source” in the Fig. 1(b)), entering the expres¬ 
sion of the Purcell factor Eq. (1) is a well-defined quantity 
in optics as well as in classical physics and antenna en¬ 
gineering. Developing this concept, in Section II we pro¬ 
pose a new methodology of calculation and measurement 
of the electric and magnetic Purcell factors through the 
input impedance of an equivalent small antenna, general¬ 
izing in Subsection IIF the results of Refs. [20, 37, 41, 42] 
and suggesting in Subsection IIB an equivalent scheme 
of nanoantenna’s Purcell effect suitable for the direct cal¬ 
culation of involved impedances. As an illustration, in 
Section HE we show how our method can be applied to 
the analytical and numerical calculations of the Purcell 
factor of nanoantennas for quantum emitters. Our ap¬ 
proach naturally results in a method for calculations and 
direct measurements of Purcell factors through the input 
impedance of arbitrary electric and magnetic dipole an¬ 
tennas. Finally, we present an example of the method 
application for radio frequencies: Section III presents an 
experimental verification of the proposed method for mi¬ 
crowave electric and magnetic dipole antennas located 
above a metal mirror. 


II. RETRIEVAL OF THE PURCELL FACTOR 
THROUGH THE INPUT IMPEDANCE 


ject. Let us attribute No 1 to the dipole and No 2 to 
the object. The total electric field created by the dipole 
1 at its origin Ei(r d ) can be decomposed into two parts 
Ei(r d ) = E n (r d ) + Ei 2 (r d ), where E n (r d ) is the field 
created by the dipole 1 in the absence of the object 2 
and Ei 2 (r d ) = E s (r d ) is the field scattered by the object. 
The total power delivered by the radiating particle to the 
environment reads 

p = E rad + Pnonrad = - ^ J Re [j { (r) • E d (r)] dV, (9) 

where V is the volume of the radiating dipole and jj is 
the electric current density in that volume. It splits into 
two parts P = P11+P12, where Pn is the power radiated 
by the same dipole in the absence of object 2 (the same 
as Po,rad as above), and 


P 12 = - - Re 


Ei 2 (r d )- / jl(r )dV 
Jv 


( 10 ) 


Here we assume that the radiating dipole 1 has a suffi¬ 
ciently small volume, such that the spatial variation of 
field E12 over V can be neglected. Since the classical elec¬ 
tric dipole moment is defined through the electric current 
density ji as 


di = — [ ji (r)dV (11) 

jw J v 

for the time dependence in the form exp (ju>t), formula 
(10) can be rewritten as 

P 12 = -|lm[dt-E 12 (r d )]. (12) 


Formula (5) after substitution of (2) for Pn gives for the 
Purcell factor F: 


tr P11+P12 , . 6^0 1 r it ( M /10\ 

F= ---= 1 + 7-T-T77— Im]^ ■ Ei 2 (r d )] (13) 
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i.e. the known result Eq.(l), which is definitely applicable 
to both classical and quantum emitters (weakly coupled 
to an arbitrary object). 

The last result can be rewritten in terms of the input 
impedances and in terms of the Green function. First, 
by definition of radiation resistance we have: 

P = Pll + Pl2 = |/l| 2 Prad = |/l| 2 (P 0 ,ra d + Rli)- (14) 


Here the effective current I\ referred to the origin r d is 
related with the dipole moment as I\ = ju>d\/l\, l\ is the 
effective length of the dipole 1. The radiation resistance 
of an optically small particle with the effective length l 
reads as [29]: 


A. General Method 


P 0 ,ra d = JT-iklf, 

07r 


(15) 


Consider an arbitrary radiating electric dipole with 
the moment d in presence of an arbitrary passive ob- 


where p = v/ro/eoUi an d k = q^fEh are the wave 
impedance and the wave number of the host medium, 
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respectively. The additional (mutual) resistance R \2 = 
R eZi 2 caused by the field scattered from the radiation¬ 
enhancing object E 12 can be found separately. This is an 
interesting and relevant problem which will be studied in 
the next subsection. 

However, formula (14) rewritten as F = P/Pn = 
Rrad/Ro,rad may already serve as a practical alternative 
to the commonly used expression (1). From the general 
theory of antennas it is well known that the input resis¬ 
tance of a short dipole equals the radiation resistance, 
when the dissipative losses inside the antenna are ne¬ 
glected [29]. Thus, if our emitter is low-loss (7 » 7 di s )> 
we can write an equivalent relation for the Purcell factor 
of an arbitrary object (inhomogeneous environment) for 
a low-loss emitter (does not matter - quantum or classi¬ 
cal): 


F = 



Ro.in 


Re Z in 
Re ^0,in 


(16) 


In some situations it may be easier to measure or calcu¬ 
late the input impedances (^o,in and Z\ n ) of the emitter 
1 in the absence and presence of object 2 than to ac¬ 
curately find the scattered field. Then formula (16) al¬ 
lows the Purcell factor through the real parts of these 
impedances. This factor for a dipole emitter is deter¬ 
mined by the modification of the resistive part of its in¬ 
put impedance. Non-radiative losses in this formula are 
present in R ln since the additional resistance R 12 is not 
purely radiative. Mutual coupling effectively brings the 
losses of the object 2 into emitter 1 . 

Let us now show that formula (16) fits another known 
representation of the Purcell’s factor - through Green’s 
function [43]. The electric field produced by a dipole di 
stretched along the z-axis is related to the dyadic Green 
function of an inhomogeneous environment G(r, r^, w) as 
follows [ 21 ]: 

E(r) = ——G zz (r, r<j, w)di. (17) 

£0 £h 

In order to relate the Green function to the input 
impedance Z m of our dipole 1 we use the Brillouin 
method of induced electromotive forces (IEMF) [29]: 

Z in =l / EiW-j^r)^. ( 18 ) 

Jv 

The value in the numerator is called IEMF in radio sci¬ 
ence, and Ii in the denominator is the current though the 
central cross section of the dipole. In the short antenna 
approximation kl\ 0 the result reads as [29]: 


_ Ei(r d )Zi 

2/i 


(19) 


This expression establishes a relationship between the 
magnitude of the input impedance of a short dipole and 
the value of the total (not only scattered) electric Held 
Ei = -E 1 Z 0 at the dipole origin. Now, recall the defini¬ 
tion of the Green function Eq.(17) and using the relation 


Iih = jwdi we obtain an expression that links the Green 
function with the input impedance of the dipole 1 : 


G zz ( 0,0, w) 


4jw£ 0 ^ 
Ijk 2 “ 


Now we may rewrite (16) in form 


Rin _ ImG Z z(0, 0,u>) 
R?? ~ ImG< 0 i(0,0, U )' 


( 20 ) 


( 21 ) 


This expression is equivalent to formula (2.4) from [43]. 

Thus, it is possible to find the Purcell factor either 
using the standard techniques, such as Eqs. (1), (21), or 
using formula (16), in terms of input resistances. 


B. Equivalent Circuit for Finding the Purcell 
Factor 

Here we explain how to find the input resistance R ln 
of an optical (e.g. fluorescent) emitter in the presence 
of an optically small resonator. Such resonators, called 
nanoantennas, are used to enhance the spontaneous emis¬ 
sion of isolated quantum emitters (see e.g. in [20, 32, 36]). 
Though this treatment is targeted to quantum emitters, 
our consideration is fully classical and based on the con¬ 
cept of electromagnetically coupled oscillators. There¬ 
fore, it is relevant to illustrate our approach by equivalent 
circuits. 

For instance, we notice that in spite of physical dif¬ 
ferences between a quantum emitter and a nanoantenna 
(which is a classical resonant scatterer), both these ob¬ 
jects in the absence of tunneling effects interact purely 
electromagnetically, and their coupling is governed by 
Maxwell’s equations. Therefore, both of them can be de¬ 
scribed in terms of resonant RLC- circuits. An attempt 
to build such schemes was done in work [37], however 
without practical results. In the present paper we intro¬ 
duce an alternative equivalent circuit for radiating sys¬ 
tems comprising an optical emitter and a nanoantenna. 
This circuit illustrates a simple algorithm for calculat¬ 
ing the additional term R 12 entering the input resistance 
Ri n in presence of object 2. This term is called mutual 
resistance R m = R\ 2 - 

First, we recall the well-known circuit model of an opti¬ 
cally small dipole scatterer excited by an external electric 
field E = z 0 E (see e.g. in [44]). The current, induced 
in a short dipole antenna of effective length l reads as 
/ = El/Z , where Z = i? rac j + Rdis + jX is the total 
impedance of the particle, see Fig. 2(a). Here we have 
split R e(Z) onto the radiation resistance i? r ad, and dis¬ 
sipation resistance i?dis- Since the induced dipole mo¬ 
ment equals di n( j = 11/joj (assume for simplicity that 
dind = z odind that holds for a spherical particle for any 
polarization and for an ellipsoidal one polarized along 
one of its axes), the inverse polarizability a -1 = E/d ln a 
reads as 

— = 77 [jkd(i? ra d + i?dis) — uiX]. (22) 

a 1 
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(a) 


(b) 


R 


rad 


R 


dis 



X= ©L-1/coC 



Figure 2. (a) An equivalent scheme of a resonant dipole scatterer. (b) Equivalent schemes of an emitter and a nanoantenna in 
terms of induced electromotive forces. 


Substituting (15) into (22), we find 

1 . k 3 .OjRdis wl 

a ^ 6TT£o£h + ^ l 2 l 2 


(23) 


This is the well-know formula for the polarizability of a 
lossy dipole scatterer which is applicable to both quan¬ 
tum emitter and nanoantenna. However, in this paper 
we neglect the induced part of the dipole moment of the 
quantum emitter as well as the hybridization of its states. 
In the weak coupling regime di(wo) = 2d. However, the 
emission spectrum has the Lorentzian shape [21], and 
this means that we have to consider the polarization of 
the nanoantenna at any frequency w. So, we use the po¬ 
larization model (23) for the nanoantenna. Note, that 
using (23) it is easy to find the general limitations on the 
absorbing and scattering cross sections of the nanoan¬ 
tenna (see e.g. a review [45]). The equivalent circuit 
of the nanoantenna is shown on Fig. 2(a) and it con¬ 
tains an IEMF E = El loaded by a series connection 
of the antenna radiation resistance R r ac j, the dissipation 
resistance i?dis; capacitive impedance 1 /juiC and induc¬ 
tive one jujL. This series connection corresponds to the 
Lorentzian model of the scatterer’s dispersion: 


1 

a 


CKO 


'( w 0 — W 2 + jujT dis) + j 


k 3 


6n£ 0 £ h 


(24) 


Comparing (23) and (24) we can relate the equivalent pa¬ 
rameters with the corresponding parameters ao, ujq and 
Tdis of the Lorentzian model: 


RrUs — l 


2 bdis 
CKO 


L = 


c = 


CKO 


Op 

l 2 ul 


(25) 


Obviously, for the resonance frequency we have Wq = 
1/(LC). 

For the emitter we also start from the general circuit 
model Fig. 2(a). This equivalent circuit corrects and re¬ 
places an incorrect scheme suggested in [37] (Fig. 3(a)). 


However, in the approximation of weak coupling as above 
we assume that the electric dipole moment di = diZo 
is fixed at any frequency corresponding to the emission 
spectrum. Because the dipole moment is related to the 
effective current of the emitter I e = jujdi/li , the equiv¬ 
alent circuit in Fig. 2(b) comprising both emitter and 
nanoantenna is driven by a fixed current source I\ = I e . 
The replacement of the circuit driven by the EMF E in 
Fig. 2(a) by the circuit shown in Fig. 2(b) which is driven 
by the current generator is granted by the well-known 
equivalent generator theorem. In Fig. 2(b) we neglect the 
dissipation in the quantum source since the main mech¬ 
anism of the decay rate is radiative (i?dis ■C -Rrad)- Here, 
for simplicity of notations, i? ra d denotes the proper radi¬ 
ation resistance of the emitter denoted above as i?o,rad- 


The IEMF describing the mutual coupling of nanoob¬ 
jects 1 and 2 in Fig. 2(b) can be replaced by mutual 
impedance Z m , whose real part comprises an additional 
radiation resistance arising in the emitter and responsible 
for the Purcell factor. The corresponding modification of 
the equivalent scheme from the mutually induced EMF to 
the mutual impedance can be accompanied by following 
speculations. The emitter induces the IEMF £ 2 1 = E 2 \l 2 
in nanoantenna 2, where the field E 2 \ is that produced 
by the emitter at the center of the nanoantenna r 2 . This 
field can be written in form E 21 = ZoA ee di, where A ee 
is the electric field of a unit electric dipole with the ori¬ 
gin at ri = evaluated at r 2 . In the case of symmet¬ 
ric mutual location of objects 1 and 2 the value A ee is 
scalar. This IEMF is related with the current induced 
in the nanoantenna as I 2 = ^ 21 /^ 2 , where Z 2 is the 
impedance of the nanoantenna. The dipole moment of 
the latter d 2 = hh/jw = E 2 R 2 .IjuZ 2 generates the scat¬ 
tered field E 12 and the IEMF E\ 2 = E 12 I 1 arises in the 
quantum emitter. Due to the reciprocity we may express 
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Mutual impedance 
1: Quantum emitter Z m =R m +jX m 

2: Nanoantenna 


(a) 





Figure 3. (a) Equivalent schemes of an optical emitter and a nanoantenna in terms of mutual impedance, (b) An equivalent 
scheme of an emitter with the mutual impedance added by the nanoantenna. 


E 12 through the same coefficient A ee : 


rewrite formula (27) as 


El2 — ^oA ee d 2 


A 2 e d\l\l2 

ju>Z 2 


(26) 


_ j^Lc If _yy-2 

1 ~ (^) “ 


(28) 


Since the current in the emitter is fixed, I\ = judi/li = 
I e , the IEMF £12 = Ei 2 l\ is equivalent to the mutual 
impedance Z m = —£i 2 /h in accordance to the equiva¬ 
lent generator theorem. The minus sign in the relation 
Z m = S 12 /I 1 appears because the IEMF £12 is directed 
oppositely to the driving current I\ (Fig. 2(b)). Thus, 
in our final equivalent scheme Fig. 3(a) the IEMF £12 
is replaced by the mutual impedance Z m describing the 
contribution of the nanoantenna into the the emitter cir¬ 
cuit. From the equivalent generator theorem and Eq. (26) 
we obtain 


Zm. — 


£12 = iVMk 

h uj 2 Z 2 ' 


(27) 


The final equivalent circuit of the radiating system 
where the presence of the nanoantenna is fully described 
by the mutual impedance Z m is depicted in Fig. 3(b). 
The radiation of the whole system is created by the cur¬ 
rent generator I\ = I e = jud/li loaded by the series con¬ 
nection of the proper impedance R la .d+ jX of the emitter 
and the mutual impedance Z m . In the reactance X of the 
emitter its proper L- and C-parameters are connected in 
series. In the mutual impedance Z m the effective mu¬ 
tual inductance L m and capacitance C m are connected 
in parallel. This difference needs to be explained. 

The input impedance of the nanoantenna is a series 
connection of resistance R 2l inductance L 2l and capac¬ 
itance C 2 - Values R 2 , L 2 and C 2 can be found from 
the Lorentzian model of the nanoantenna - formulas 
(25). Substituting Z 2 = R 2 + ju)L 2 + l/jwCz, denot¬ 
ing uj 0 = I/VC 2 L 2 , and assuming that ui « ujq, we may 


It is the standard formula of the circuit theory describ¬ 
ing the impedance of a voltage transformer loaded by 
a low-loss parallel circuit resonating at too- In this 
formula L e g = C 2 ^o/so is the effective inductance of 
the parallel circuit and the dimensionless value N = 
y/£olil 2 A ee /u)^/JIo is an effective transformer parameter 
(called turns’ ratio in the electrical engineering). In the 
vicinity of the resonance the dispersion of Z m is mainly 
determined by the denominator and we may neglect 
the frequency dependence of the effective transformer 
putting in Eq. (28) N rts y/£olil 2 A ee /u}o^/JIo. Then for¬ 
mula (28) describes the impedance of a parallel circuit 
with mutual inductance L m = fioC 2 N 2 / Eq and mutual 
capacitance C m = £oL 2 /N 2 fi 0 connected to effective re¬ 
sistors (nonzero ones in both capacitive and inductive 
branches) which are responsible for the mutual resistance 

Rm- 

The value R m = R 12 - the real part of the right-hand 
side of Eq. (28) comprises both radiative and dissipative 
resistance added to that of the emitter due to the pres¬ 
ence of a nanoantenna. In the quasi-static approximation 
the value of A ee is real. Then N is real and positive that 
results in the Purcell effect larger than unity. If N 1 
the Purcell factor at the resonance frequency may take 
huge values. 

Since the driving current is fixed, the power deliv¬ 
ered by the emitter to its environment is equal to P = 
|/i | 2 (7? r ad + Rm)• The Purcell factor in accordance to 
(16) takes the form 


F = 1 


Re Z r 
R 


= 1 


rad 


Gnl^l 2 

rjui 2 k 2 l\ 


Re 


A 2 

) , (29) 
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where we have used formula (15) for i? ra d and sub¬ 
stituted relation (27). Now, applying the model of a 
Lorentzian scatterer to the nanoantenna we may express 
the impedance Z 2 of the nanoantenna through its po¬ 
larizability a 2 = oina- Really, g ?2 = E^l^/juiZ 2 and 
a 2 = d 2 /E 2 \. Therefore, (29) can be rewritten as 

f\'jr 

F = I + -3 - Re(ja 2 A 2 ee ). (30) 

uj 3 rj£h v ee ' 

This expression clearly shows that the Purcell factor does 
not depend on the emitter 1 - only on the nanoantenna 
2 and their mutual location. Therefore we speak on the 
Purcell factor of an object at a point with radius vector 
iq — r 2 with respect to the object. This factor is applica¬ 
ble to an arbitrary dipole emitter located at this point. 

Now, it is time to make some other important com¬ 
ments. First, the problem of mutual coupling which 
we have solved above corresponds to the steady regime 
and is self-consistent at every frequency. The Purcell 
factor has the physical meaning at frequencies close to 
wo, since the emission has a finite decay rate and its 
spectrum has nonzero bandwidth. Second, our analy¬ 
sis keeps valid in the case when the nanoantenna 2 has 
the resonance frequency w 02 different from the emission 
one wo = W 01 . Still formula (30) holds and the equiva¬ 
lent scheme depicted in Fig. 3 remains adequate, but the 
mutual impedance Z m is not anymore that of a simple 
parallel circuit connected through the transformer. How¬ 
ever, if the difference between wo and w 2 0 is large, Ima 2 
becomes too small at the emission frequency wo, and the 
Purcell factor is close to unity. 

The second comment is more important. In fact, the 
factor A ee (electric field of a unit dipole with origin iq 
evaluated at r 2 ) is complex due to the retardation effect. 
Its imaginary part is relevant for calculation of the Pur¬ 
cell factor at the frequencies different from the resonance 
frequency of object 2. Moreover, it is not exactly deter¬ 
mined by the field of a unit dipole at the geometric center 
of the nanoantenna r 2g . The electromotive force induced 
by a point emitter in nanoantenna 2 may be found ac¬ 
curately - via the integration of the local field E 2 i(r) 
over the volume of the nanoantenna. If the local field 
is strongly non-symmetric with respect to its geometric 
center, the effective center r 2 of the nanoantenna shifts 
from the point r 2g towards the emitter. Further, from the 
classical antenna theory [29] it is known that for a two- 
element array of dipole antennas the mutual resistance is 
positive only when the antennas are collinear. This mu¬ 
tual location of dipoles 1 and 2 corresponds to Fig. 4(a) 
when the dipole moment of the emitter is stretched ra¬ 
dially towards a plasmonic nanosphere. In this case for 
small distances G between the emitter and the sphere we 
may approximate lmA ee = 0 and A ee « 1/27T£o £hD 3 , 
where D = a + G. In accordance to Eq. (30) this results 
in the Purcell factor higher than unity. However, if the 
dipole is located with respect to the nanosphere so that 
their dipole moments are parallel and not shifted, the 
interaction of dipoles becomes destructive. In this case 


one cannot neglect ImA ee , moreover, its contribution for 
distances D comparable with l\ and l 2 significantly ex¬ 
ceeds that of the real part [29]. Then the second term 
in (30) becomes negative and makes the Purcell factor 
smaller than unity. This corresponds to the known situ¬ 
ation: the mutual resistance R m of a transmitting dipole 
and a closely located reflector antenna is negative, and 
the enhancement of the directionality is accompanied by 
the decrease of the efficiency [29]. For this case the an¬ 
tenna theory gives i? ra d < |i? m | [29]. The energy balance 
is respected and we have F > 0. 

Last comment of this subsection refers to the approxi¬ 
mation of the fixed dipole moment di = d. This approxi¬ 
mation restricts the allowed strength of the dipole-dipole 
interaction by relatively modest values of the Purcell fac¬ 
tor. Very high values of this factor would imply a very 
strong interaction of quantum object 1 (with both elec¬ 
tromagnetic field and with nanoantenna 2). For this very 
strong coupling the approximation of predefined eigen¬ 
states becomes inadequate. Then, one has to solve the 
self-consistent problem for eigenmodes of the radiating 
system comprising the emitter and the nanoantenna. The 
solution within the framework of the semi-classical the¬ 
ory results in the so-called Rabi splitting [46]. Instead 
of one emission frequency wo two emission frequencies 
arise, corresponding to two (one in-phase and one out- 
of-phase) spontaneous oscillations in the resonating and 
radiating system formed by objects 1 and 2. This ef¬ 
fect has a well-known analogue in the classical theory of 
two coupled oscillators and can be of course described by 
a corresponding equivalent circuit. However, this phe¬ 
nomenon is beyond the framework of the present paper. 
Similarly, we do not consider the case of strong coupling 
when W 02 7 ^ woi- 

C. Validation of the Equivalent Circuit 

To validate our circuit model we apply it to an ex¬ 
plicit structure depicted in Fig. 4(a). First, let us show 
that formula (30) based on the equivalent circuit fits the 
known analytical solution [47]. In that paper the Purcell 
factor was calculated using the exact solution of the elec¬ 
trodynamic problem of a dipole radiating in the presence 
of a sphere of arbitrary radius a filled by an isotropic 
material of (generally complex) permittivity e s . Formula 
( 6 ) of that paper refers to the radial polarization of the 
dipole and its location outside the sphere. It is a series 
in which the first term corresponds to the dipole polar¬ 
ization of the sphere, i.e. for optically small spheres we 
may neglect the other terms. 

In the framework of this approximation formula for the 
radiative Purcell factor [47] reads: 

F w 9\ji(kD) + 6 1 /4 2 ) (H?)|7(/cD) 2 , (31) 

where j\ (X) and h\ (X) are, respectively, spherical 
Bessel’s and Hankel’s functions with n = 1, D = a + G 



(a) 


(b) 



1: Fluorescent 
quantum dot 


2: Au nanosphere 

plasmon resonance 



2: Si nanosphere 

Magnetic Mie resonance 



Figure 4. (a) A fluorescent emitter over a plasmonic (e.g. golden) nanosphere has F > 1 when its dipole moment is radially 
directed since in this case A ee > 0. (b) The same emitter over a dielectric (e.g. silicon) nanosphere has F > 1 when its dipole 
moment is azimuthal: in this case A ern > 0. 


is the distance between the emitter and sphere centers, 
k = q^feh is the wave number of the host medium and 
coefficient b\ is given by formula [47]: 

b _ £hji(ka)[k s aj 1 (k s a)]' - £ s j 1 (k s a)[kaj 1 (ka)} 1 

e s ji(k s a)[kah^\ka)]' - Ehh i {\ka)[k s aji(k s a)}' 

(32) 

Here k s = q^fel is the wave number inside the sphere. 
The dipole approximation is valid when |fc s |a -C ^ r, prac¬ 
tically when |fc s |a < 1. To compare the radiative Purcell 
factor (31) with our result (30) we have to remove losses 
that automatically equates the total Purcell factor to the 
radiative one. Therefore, we assume that e s is real. Then 
k s is either real (if e s > 0) or imaginary (if e s < 0). In 
both these cases the following approximations are suit¬ 
able for the spherical functions entering (32): 


X 

, X = k s a,X = ka. 

(33) 

With these substitutions the differentiation in (32) be¬ 
comes elementary, and we obtain for b\ the result bi = 
jB , where B is a real value: 


B 


2 (ka) 3 £ a - e h 
3 £ s + 2 Eh 


(34) 


Let us restrict the analysis by the case kD -C it . Then, 
the asymptotic relations (33) are suitable for X = kD. 
Formula (31) with substituted expressions (33) and (34) 
can be rewritten in a form 


F 


|1 +3 


3 b\ 
(kD) 3 


= 1 + 


9 B 2 

q 6 D 6 £ 3 ' 


(35) 


Substitution of (34) into (35) results in 


Our circuit model resulted in formula (30) which can 
be rewritten as 

F = 1 - 3 n< L- Im ( a 2 A 2 ee ) . (37) 

U> 3 y/£h 

The quasi-static approximation for A ee has been already 
introduced: 


2TTEoE h D 3 ' 


(38) 


Since the sphere is lossless, we have in accordance to (24): 


Im a .2 





(39) 


The quasi-static polarizability «Qg of a small sphere is 
well-known (see e.g. in [48]), and we have: 


«2 ~ «qs = 4na 3 E 0 £h £ ' b - " h ■ (40) 

e s + 2 £ h 

Substituting (38), (39) and (40) into (37) we obtain for¬ 
mula (36). Thus, the strict electrodynamic model and 
the present circuit model meet one another within the 
framework of the dipole approximation. 

In order to validate our circuit model for a more in¬ 
teresting case when the sphere is resonant (plasmonic 
nanoantenna) we consider an explicit structure of a 
golden sphere of diameter 2a = 40 ran. The radially po¬ 
larized emitter is located at the distance G = 10 nrn from 
its surface. Values of the permittivity are taken from the 
experiments of Johnson and Christy [49, 50]. The radi¬ 
ating system is located in the air £h = 1. In Fig. 5 we 
present our calculation of the Purcell factor performed 
with the use of (37) in comparison with the Mie theory. 
In our calculation we have complemented formula (40) 
by radiation losses in accordance to (23): 


F 


1 


4a 6 / e s - Eh \ 
D 6 \£ s + 2 £ h ) 


(36) 


Oil 


1 k 3 \ 

OQS '^67 T£ 0 £ h ) 


-1 
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For the plasmon resonance band our model is in agree¬ 
ment with the exact calculation. Our rough approxima¬ 
tion for A ee works in this band because at the resonance 
the dipole eigenmode is realized. The excitation mecha¬ 
nism is not very important, and the sphere is polarized 
by an emitter as if it were excited by a plane wave - 
nearly uniformly. The model becomes less accurate be¬ 
yond the resonant band, where strong non-uniformity of 
the external field E 21 implies strong non-uniformity of 
the polarization decaying versus the distance from the 
emitter. Due to this decay the origin r 2 of the dipole 
d 2 shifts towards the emitter, the effective distance de¬ 
creases compared to D and A ee increases compared to 
(38). Therefore, it is not surprising that our model uti¬ 
lizing the simple approximation (38) underestimates the 
Purcell effect at low frequencies. 


D. Extension of the Circuit Model 

Next, let us extend the equivalent scheme and general¬ 
ize formula (30). First, let us see that the equivalent cir¬ 
cuit keeps valid if the nanoantenna is a magnetic dipole. 
Of course, we mean artificial magnetism when the vortex 
polarization currents in the subwavelength particle result 
in its magnetic dipole moment. Qualitatively, this insight 
is applicable, for example, to a submicron silicon sphere 
at its magnetic Mie resonance (see e.g. in [51-53]). In 
Fig. 4(b) we have depicted the corresponding radiating 
system. The z-directed electric dipole d\ of the emitter 
1 induces at the center of the nanosphere 2 a magnetic 
dipole m 2 = xo?7i2 which is related to the local magnetic 
field via the magnetic polarizability /?2 = m. 2 /H x . In this 
definition H x is the local field acting on the magnetic 
dipole. Here H x = H 21 is the magnetic held produced by 
the electric dipole di in the plane 2 = 0 at the distance 



Wavelength, nm 


Figure 5. The radial Purcell factor of a golden nanosphere 
of diameter 40 nm at the distance G = 10 nm in air: our 
circuit model (red solid curve) and exact Mie theory (blue 
dashed curve). Values of the permittivity are taken from the 
experiments of Johnson and Christy [49, 50]. 


D. It can be written as H 21 = jujA ern d\ , where in the 
quasi-static limit A em « l/47rZ? 2 . 

The magnetic dipole antenna can be modeled as an 
optically small loop with an effective area S and effec¬ 
tive electric loop current I which is considered uniform 
around the loop. The magnetic dipole moment is equal 
to m = /io Sin, where n is a unit vector to the loop 
plane. The input impedance of the effective loop antenna 
equals to the ratio of IEMF £ = jui^H n S (where H n is 
the normal component of the local magnetic field) to the 
electric loop current I. So, in the present case the IEMF 
for the magnetic nanoantenna 2 resulting in the mag¬ 
netic moment m 2 is equal to £21 = jw/zoL^iiS, where S 
is the effective area of the polarization current loop of the 
nanosphere (it will cancel out in the result). The induced 
magnetic moment m 2 = Hq£\ 2 S/Z 2 comprises the factor 
w 2 in the magnetic analogue of (22): 

IT = Tq 2 - uX 2 ). (41) 

P2 uPS 

Respectively, the Lorentzian model of the magnetic po¬ 
larizability of a scatterer differs from the model of the 
electric polarizability by the factor ui 2 (see e.g. in [51]). 
The analogue of expression (24) takes form: 


1 

0 


1 

w 2 /3 0 


Wn 


w 2 + jcuTdis) + j 


fc 3 

67rpo' 


(42) 


All other formulas of the Lorentzian model keep valid. 

Accordingly, the equivalent circuit remains applica¬ 
ble. Magnetic moment m 2 produces the electric field 
E 12 = juA em m 2 , comprising due to the reciprocity the 
same coefficient A em which enters H%i. The correspond¬ 
ing IEMF £12 = E 12 I 1 is recalculated into the mutual 
impedance in the same way as above. Reproducing the 
same steps as for the electric dipole nanoantenna we come 
to the parallel-circuit formula (28) for Z m with substitu¬ 
tion N = LuS 2 hA em /c for the transformer parameter. 
For the Purcell factor we obtain an analogue of (30) in 
the form: 


F=l + — Re(jP 2 A 2 ern ). (43) 


If the response of the nanoantenna comprises both elec¬ 
tric and magnetic dipoles, each of these modes is de¬ 
scribed by its own equivalent scheme. If these dipole 
moments resonate at the same frequency both equivalent 
circuits are similar and can be unified. A more com¬ 
plicated equivalent scheme would correspond to different 
resonances of the electric and magnetic modes. However, 
it is important to notice that both the electric and mag¬ 
netic modes obviously contribute into the total mutual 
impedance, and both these contributions can be con¬ 
structive. So, the excitation of an additional mode in 
the nanoantenna may increase R ra d enhancing F. The 
same refers to higher multipoles of the nanoantenna: each 
of the multipole modes contributes into total Z m , and 
the coinciding or closely located resonances of high-order 
multipoles may result in huge values of the Purcell factor. 
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Wavelength, nm Wavelength, nm 


Figure 6. The Purcell factor extraction through a change of the input impedance in optics, (a) Illustration of a point dipole 
source located close to the dielectric spherical nanoparticle of the radius a = 70 nm. (b) and (c) Purcell factor dependence on 
the emission wavelength for the parallel (b) and perpendicular (c) dipole orientation with respect to the sphere. 


E. Theoretical Verification of the General 
Approach 

First, let us notice that our general approach resulted 
in formula (16) is a useful alternative to the conventional 
methods of calculating the Purcell factor. Although vari¬ 
ous numerical methods to solve the problems of nanopho¬ 
tonics and metamaterials have become widespread [54], 
direct numerical calculation of the Purcell factor using 
Green’s function technique (21) or scattered field tech¬ 
nique (1) faces fundamental difficulties. Indeed, the ex¬ 
act calculation of the microscopic field inside the quan¬ 
tum emitter as well as the exact calculation of the Green 
function at this point is challenging and time-consuming. 
Next, as shown in Ref. [55] another known method of the 
Purcell factor calculation through the volume and qual¬ 
ity factor of the cavity mode (see e.g. in works [4, 10, 56]) 
gives a strong disagreement with the accurate theoretical 
model (1), especially for plasmonic nanostructures. In fi¬ 
nite systems and systems without losses the method of 
integrating the radiated power flow through some spheri¬ 
cal surface surrounding the radiating system has become 
popular. However, in structures with losses this method 
depends on the choice of the integrating sphere (even low 
losses may strongly deviate the result since the integra¬ 
tion surface is very large). Finally, all these methods can 
not be realized experimentally and extended to the ra¬ 
dio frequency range (which is one of the purposes of the 
present study). 

In commercial software packages, such as CST Stu¬ 
dio [57], a point dipole can be modeled as an optically 
very short dipole of a perfectly conducting wire excited 
by an ideal current source. Since it has a finite length 1 1 , 
this dipole 1 in free space has a certain finite impedance, 
whose real part I?o,in is its radiation resistance. In pres¬ 
ence of an arbitrary object 2 the IEMF £12 arises in the 
dipole and its input resistance modifies i?; n ^ i?o,in- The 
input impedance results from exact simulations with the 
use of any reliable commercial software. The result for 
the Purcell factor F = Rin/Ro,in should not depend on 
the length l\ of the equivalent Hertzian dipole. This 
method appears to be very practical and convenient for 


nanooptics. Moreover, it is more universal than all the 
aforementioned methods, because it is equally applicable 
to systems with or without losses. 

To validate the general formula (16) we have studied 
the structure depicted in Fig. 4(b). In Fig. 6(a) the ge¬ 
ometry of the problem under consideration is recalled. 
The quantum source is modeled in CST Microwave Stu¬ 
dio as a Hertzian dipole of length 10 nm. The dielectric 
spherical nanoparticle of radius a =70 nm and relative 
permittivity 15 is located at distance G from the dipole. 
The Purcell factor retrieved from numerical simulations 
as F = Rin/Ro.in is compared with the exact solution [47] 
in which now the series has been accurately evaluated. 
We studied both parallel and orthogonal dipole orienta¬ 
tions, corresponding to Figs. 6(b) and (c), respectively, 
for three values of G. The exact solution and our results 
are in excellent agreement. For the orthogonal orienta¬ 
tion at wavelength A ~ 570 nm the sphere experiences the 
magnetic Mie resonance, and at A « 390 nm - the electric 
dipole and magnetic quadrupole (makes the largest con¬ 
tribution) Mie resonances. Unfortunately, our simplistic 
model resulting in formulas (30) and (43) does not of¬ 
fer enough numerical accuracy due to two factors. First, 
the electric dipole mode cannot be neglected at the mag¬ 
netic resonance. Second, the electric resonance holds at 
higher frequency, where the electromagnetic response of 
the nanosphere is not purely dipolar. However, for our 
current purpose it is enough that the exact version of our 
method ■ formula (16) - gives an excellent accuracy. 


F. Purcell Factor for Radio Antennas 

Now, let us go beyond the optical frequency range and 
extend the whole concept to radio frequencies, includ¬ 
ing microwaves, millimeter waves and teraherz frequency 
range. Instead of a quantum emitter let us consider a 
dipole antenna 1 interacting with an arbitrary object 2, 
as it is sketched in Fig. 7. If the dipole 1 is resonant, e.g. 
has length l\ = A/2 it can be excited by a short pulse (an 
analogue of the optical pumping) and will irradiate its en¬ 
ergy at its resonant frequency ujq during the finite emis- 
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Figure 7. A schematic illustration of the radiative dipole an¬ 
tenna placed near an arbitrary scattering object. Each cur¬ 
rent density element of the antenna interacts with itself and 
other elements of the current (En), as well as with an object 
(E 12 ). 


sion time I/70. If the radiation quality of the antenna 
is high the time I/70 is very long in terms of the period 
27r/wo- It may be reduced to I/7 <C I/70 if an object 2 
is located in the vicinity of the antenna 1 which increases 
its radiation resistance. Object 2 is not obviously a res¬ 
onator tuned to the same frequency as it is adopted in 
optical applications of the Purcell effect. In accordance 
to our consideration in subsection it can be an arbitrary 
object constructively interacting with the antenna. Then 
the general equivalent circuit shown in Fig. 2(b) remains 
valid and results in the mutual impedance Z m . Of course, 
in the general case the mutual impedance is not obviously 
that of a parallel RLC- circuit. What is essential that 
R e Z m = Ri 2 should be positive and increase the input 
resistance of antenna 1. The impressed current I e = Ii is 
then determined by the dipole moment of the antenna 1 
at the moment when the external pulse ends. In fact, this 
consideration, and the representation of the object 2 via 
the mutual impedance Z m have been known in antenna 
engineering for a long time, see e.g. [58]. 

It is difficult to significantly increase the radiation re¬ 
sistance of an already efficient antenna - that with the 
resonant length 1 1 = A/2. Absolute values of R m may be 
noticeable in this case, but the relative contribution will 
be modest. The concept of the Purcell factor becomes 
relevant for a short dipole - that with a low radiation re¬ 
sistance Ro.radi much lower than the internal resistance 
of the voltage generator applied to the radio antenna. As 
a rule, this is the output resistance of the feeding trans¬ 
mission line which usually equals R out = 50 Ohms. If 
-Ro.rad "C -Rout, the presence of a low-loss object insert¬ 
ing positive mutual resistance may lead to much better 
matching of the effective generator to the antenna and 
therefore to much higher radiation. At first glance, this 
radiation gain has nothing to do with the Purcell fac¬ 


tor, which describes the emission regime. However, for a 
very short dipole l\ <C A/2 these values are equal to one 
another. 

The proper reactance of a short dipole antenna is 
capacitive. The spontaneous emission (quasi-harmonic 
radiation after a short pulse) is possible if the output 
impedance of the feeding line has the inductive reactance 
connected in series with 7? out . Then, in the absence of 
object 2 the emission is still described by the current 
source R, at the frequency ujq = 1/y/LC loaded by the 
resistance of the feeding line R out = 50 Ohms and the 
radiation resistance Ro, ra d <C R ou t- Most part of the en¬ 
ergy is lost in Rout and only a small portion of the pulse 
energy is irradiated. The presence of object 2 changes 
this distribution increasing the radiation resistance and 
the decay rate multiplies by F = R ra d/Ro,rad- 

In the regime of the usual transmission at the fre¬ 
quency too, the steady-state voltage V at the output 
of the feeding line is loaded by the resistance R out = 
50 Ohms and the antenna input impedance Z- m . The 
input impedance of a small antenna consists of a small 
radiation resistance Ro, ra d Rout and a very high re¬ 
actance X. In this case the current I\ = V/(R out + 
Ro,rad + jX) « V/{R ou t + jX) <C V /Ro, ra d and only 
a small portion of the supplied power is radiated. The 
power is mainly reflected from the antenna back to the 
generator. The presence of object 2 increases R ra d, he. 
improves the matching of the antenna to the feeding 
line. The radiated power increases in accordance to 
formula P ra( j = |/i| 2 R ra d- However, matching remains 
poor since R ra d ■C |R 0 ut + jX\. Therefore, we can write 
h = V/(Rout + R ra d +jX) « V/(R out + jX). The radi¬ 
ating current does not change in the presence of object 
2 though the input resistance of the antenna 1 changes! 
The increase of the radiated power is solely described by 
the increase of the radiation resistance. Therefore, the 
gain in the transmitted radiation is equal to the Purcell 
factor F — Rrad/Ro.rad* 

Briefly, for a very poor transmitting antenna 1 we may 
find the Purcell factor of object 2 from the usual radiation 
gain of the same antenna 1 in presence of the object 2. 
This factor describes the emission of the pulse energy by 
antenna 1 in presence of the radiation-enhancing object. 
It does not depend on the antenna itself and is fully de¬ 
termined by the properties of the object and its location. 
Vice versa, we can predict how much the antenna will ra¬ 
diate due to the presence of the object if this antenna is 
tuned into resonance, excited in the absence and presence 
of the object by a pulse voltage, and find the decay rate 
of its emission after the pulse is gone. We should stress 
that the Purcell factor of object 2 measured with the use 
of an antenna is not the same as the radiation enhance¬ 
ment of this antenna in the presence of object 2. Only 
in an important special case when the probe antenna is 
a very poor emitter they are approximately equal. The 
observation of this equivalence dramatically extends the 
field where the notion of the Purcell factor is relevant. 

Our last extension concerns the Purcell factor of an ar- 
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bitrary object acting on a magnetic dipole antenna. We 
have already noticed that the magnetic dipole antenna 
is an optically small loop (can be multi-turn [29]) with 
an effective area S and electric current I which is prac¬ 
tically uniform around the loop. The magnetic dipole 
moment m = /xoSTn, where n is a unit vector to the 
loop plane is related to the effective magnetic current as 
/ TO = to = —juim. The input impedance of the loop an¬ 
tenna equals to the ratio of IEMF £ = jojfj, 0 H n S (where 
H n is the normal component of the local magnetic field) 
to the induced electric current I and can be rewritten as 
^in,m = Hn/Im- This offers a full analogy with the elec¬ 
tric dipole antenna and corresponds to the duality princi¬ 
ple. It is clear that the input impedance of the magnetic 
antenna is related to the Green function at the magnetic 
dipole origin as G zz (0, 0,w) = q 2 . After ex¬ 

tracting the imaginary part from the last expression, we 
obtain the Purcell factor in the form (16) [41]. So, all the 
theory developed above including the equivalent circuits 
keeps valid. 


III. MEASUREMENT OF THE PURCELL 
FACTOR IN THE MICROWAVE SPECTRAL 
RANGE 


Now let us demonstrate the application of our method 
experimentally, retrieving the Purcell factor from mea¬ 
sured input resistance of a radio antenna using Eq. (16). 
The input impedance of an antenna can be easily deter¬ 
mined from the S'-parameters. Namely, for a dipole an¬ 
tenna connected to a one-mode waveguide (e.g. a coaxial 
cable), the quantity R\ n is related to the reflection co¬ 
efficient Sir measured at the waveguide input and the 
characteristic impedance of the waveguide [29]: 


1 - [Re(Sn)] 2 - [Im(Su)] 2 
[1 - Re(S'n)] 2 + [Im(iSii)] 2 ' 


(44) 


In our experimental verification of the technique ob¬ 
ject 2 is a flat copper plate of optically large size and 
the antenna 1 is located near its center. This plate in 
the microwave range emulates the perfectly conducting 
plane and the Purcell effect in this case may be referred 
to as a special case of spontaneous emission near an in¬ 
terface [59-70]. For the perfectly conducting interface 
a simple analytical result for the Purcell factor was ob¬ 
tained in [68, 71]. The expression for the electric ( F e ) 
and magnetic ( F m ) Purcell factor for either parallel (||) 
or perpendicular (_L) orientations are as follows: 
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(45) 


where r) = 2 qh, h is the height of the (electric or mag¬ 
netic) dipole above the metal. The upper sign corre¬ 
sponds to an electric dipole, the lower to the magnetic 


one. We have compared the predictions of Eqs. (45) with 
F = it’jn/it’o.in- The experimental setup is schematically 
shown in Fig. 8. As two stems of an electric dipole an¬ 
tenna we use brass wires of length 0.4 cm soldered to the 
internal and external veins of the coaxial cable connected 
to a vector network analyzer. The wave impedance of the 
cable is equal to Z w = 50 that guarantees the regime 
Ro.rad Z w . Magnetic dipole source is realized as a 
wire ring with the diameter 1 cm connected similarly. 
The measurement is performed in the spectral range 5- 
14 GHz which corresponds to wavelengths from 2.14 to 
6 cm. The object 2 was a polished stainless steel sheet 
with sides 180 x 210 cm (the smallest mirror side greatly 
exceeds the largest wavelength and the diffraction effects 
are negligible). The antennas were attached to an arm of 
a precise coordinate scanner which moved in the vertical 
directions, allowing us to measure the Purcell factor as 
a function of the emitter height. The main experimental 
results for electric and magnetic antennas are shown in 
Fig. 8b,c (squares and triangles correspond to two orien¬ 
tations of the electric and magnetic antennas). The solid 
blue and dashed red curves represent the theoretical val¬ 
ues of the Purcell factor (45). Experimental and theoret¬ 
ical results are in excellent agreement. The Purcell factor 
exhibits oscillations with the period on the order of the 
wavelength when the source is moved vertically. These 
oscillations are due to the interference pattern which ex¬ 
hibits in the radiation resistance R[ n « i? ra d- It clearly 
indicates that our general formula (16) is applicable far 
beyond the quasi-static interaction between objects 1 and 
2 assumed in the previous section. When h increases F 
eventually saturates at unity. 

Slight disagreement can be noticed for the magnetic 
antenna. It is explained by a slight current inhomogene¬ 
ity around the ring. This inhomogeneity appears when 
the magnetic dipole is parallel to the metal plane i.e. the 
loop is in the vertical plane. Definitely the lower half 
of the loop is stronger capacitively coupled to the metal 
plane than the upper one, and it results in this inho¬ 
mogeneity. Notice, that for very small h we could not 
measure the Purcell factor due to the finite size of the 
our antennas when Eqs. (45) become inapplicable. Inter¬ 
esting, that the deviations from these formulas were also 
observed for a quantum dot, located too closely to the 
mirror [24]. 

Importantly, in the microwave frequency range the 
electric dipole antenna is usually fed by a coaxial cable 
whose thickness is not negligible, though optically very 
small. This factor results in the radiation from the cable 
open end and affects the measured Purcell factor. We 
directly measure not the input resistance i?i n of the an¬ 
tenna but the sum of Ri n and SR, where the last term 
in the radiation resistance of the open cable. Therefore, 
we have separately measured the input resistance of the 
open end of the cable which obviously equals to SR and 
subtracted it from R ln found with the use of Eq.(44). 
Otherwise, the disagreement in Fig. 8b,c would be more 
noticeable. 












13 



electric dipole 

(a) 


electric dipole magnetic dipole 




h/X 


h /x 


Figure 8. (a) Geometry of the experiment to measure the value of the Purcell factor for an electric dipole antenna near a 
perfect metallic mirror, (b) Measured results for the Purcell factor (symbols) along with the analytical results Eq. (45) for 
parallel and perpendicular orientations of the electric dipole antenna with respect to the mirror, (c) The same results for the 
magnetic dipole antenna. 


CONCLUSIONS 


In this paper, we have analyzed theoretically and ex¬ 
perimentally the classical counterpart of the Purcell effect 
for subwavelength electric and magnetic dipole antennas. 
We have generalized the approach accepted in nanopho¬ 
tonics to the case of microwave antennas and recovered 
the known expression for the Purcell factor via the imagi¬ 
nary part of the electromagnetic Green’s function. Using 
this result, we propose a new method to directly measure 
the Purcell factor through the input impedance of small 
antenna. We have experimentally verified the technique 
for both electric and magnetic dipole antennas. The tech¬ 
nique has been also successfully applied to a cornerstone 
problem of the all-dielectric nanophotonics: Purcell effect 
due to the Mie resonances of dielectric sphere. We be¬ 


lieve that the proposed Purcell factor extraction method 
is versatile and can be used in various frequency ranges: 
from radio to optics. 
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